Inelastic Collapse of Three Particles 
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A system of three particles undergoing inelastic collisions in arbitrary spatial dimensions is studied 
with the aim of establishing the domain of "inelastic collapse" — an infinite number of collisions 
which take place in a finite time. Analytic and simulation results show that for a sufficiently small 
restitution coefficient, < r < 7 — 4y/3 w 0.072, collapse can occur. In one dimension, such 
a collapse is stable against small perturbations within this entire range. In higher dimensions, 
the collapse can be stable against small variations of initial conditions, within a smaller r range, 
< r < 9 - 4^ » 0.056. 

PACS: 47.50.+d, 05.20.Dd 



I. INTRODUCTION 



A system of particles interacting only through inelas- 
tic collisions is a useful idealization of granular materials, 
and has been much investigated recently |i] jlC|]. Inelas- 
ticity can make such a system evolve into a collapse state, 
in which several of the particles collide an infinite number 
of times in a finite time interval. 

Inelastic collapse in one dimension is well understood 
"~|. In two dimensions, McNamara and Young carried 
out numerical investigations and found some evidence for 
inelastic collapses of three particles ||. To understand 
the collapse mechanism in higher dimensions, we study 
the behavior of three particles in a particular model. Our 
model involves collisions which preserve the total momen- 
tum and the components of the momentum perpendicu- 
lar to the line of centers. The component of the relative 
velocity along the line of centers is reversed (as in an elas- 
tic collision) and reduced by the restitution coefficient, 
< r < 1. We look at a situation in which one particle 
(labeled zero) takes part in all collisions. The other two 
particles (labeled one and two) are alternatively a collider 
and a spectator. We assume that all particles have the 
same mass, and that particles one and two have identical 
radii. 

There are two possible reasons that previous numeri- 
cal studies might have shown collapse. One scenario is 
that the collapsed state is represented by one or many at- 
tractive fixed points, so that the collapsing orbit can be 
stable against small variations in initial data. The other 
scenario is that each orbit is unstable but that the infin- 
ity of collapsing orbits produces an observable collapse 
probability. For the specific example of three particles, 
we find an attractive fixed point for all dimensions and 
sufficiently small coefficient of restitution, r. Thus, we 



establish the possibility of the first scenario. The second 
scenario is still possible, but we have seen no evidence for 
it. For a larger r, there is an interval in which the fixed 
point is unstable against changes in the initial conditions. 

For dimensions greater than one, after the collapse has 
occurred, the particles can separate from one another. 
Thus beyond d = 1, inelastic collapse is an event, not an 
end-point, in the "lives" of the particles. 

To study the collapse, we use the methods of dynami- 
cal systems theory. Specifically, we examine the situation 
in which particles one and two are very close to particle 
zero and aimed so that the system is very close to the 
collapsing fixed point. Figure 1 shows a typical config- 
uration. Particle keeps colliding with particle I and 
particle 2 repeatedly, and an inelastic collapse may oc- 
cur. After many collisions, the distances between the 
particles become small, and the remaining collisions take 
place so rapidly that the relative motion of the particles 
is small. Therefore, in the inelastic collapse, the angle 
9 approaches a limiting value as the number of collision 
goes to infinity. We identify the constant-0 fixed point 
and then investigate its stability. We find that a fixed 
point of an inelastic collapse exists only when the final 9 
obeys: 

4 ^ (i) 



cos 6* > 



1 + r 

A collapse state will occur whenever this criterion holds 
and also steric effect do not block off the required col- 
lisions. (For example such a blockage will always occur 
at 9 = 7T.) The stability analysis implies that for the 
collapse to be stable against small perturbations in the 
initial velocities, a stronger condition is required, namely 



cos 9 > 



2^f(l+ ^r) 



(2) 
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II. THE COLLISION MODEL 



III. FLAT SURFACE APPROXIMATION 



We use the standard model of inelastic collision: due 
to a collision the component of the relative velocity of the 
colliders along the line of centers, changes by a factor of 
— r. We denote by Uj and Xj the velocity and the position 
of the jth particle at the instant before a collision occurs. 
Let us consider a collision between particles 1 and 0. In 
the course of the collision, the velocities of the particles 
change to: 



{![ = ui — A, 
u Q = u Q + A, 
u 2 = u 2 . 



(3) 



Here, the momentum transfer is A. It must point in the 
direction of the line of centers. In terms of the coefficient 
of restitution, r, this transfer is given by the expression 
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-(x 1 - x )[(xi - x ) ■ (ui - u )}. 



(4) 



Here we have assumed that the radii of the colliding par- 
ticles sum to unity so that, at the point of collision: 



(xi - x ) 2 = 1. 



(5) 



To do dynamical systems theory, we wish to look at the 
very same process repeatedly. Therefore, we introduce 
the superscripts, c, denoting the collider and s denoting 
the spectator particle, as well as a subscript, n, to denote 
the instant before the nth collision occurs (Figure 1). For 
simplicity, we take the velocity and the position of parti- 
cle to be zero. In order to make sure that the velocity 
of particle continues to vanish after the collision, we 
view the post-collision system from a frame moving with 
velocity u' = A. Then equations (||) - @ read 



U r, 
ill. 



2A, 
A. 



(6) 



These equations are supplemented by the conditions: 
A = i±Vj(^-<), 



1. 



(7) 



(8) 



Additionally, both the velocity and the position of parti- 
cle vanish. Notice that a collider becomes a spectator 
immediately following a collision. The above equations 
are complemented by the equations corresponding to the 
positions of the particles at the next collision, 



L n+1 



Xj, ~\~ tjill 



Xj, ~\~ tjill 



(9) 



The time interval between the nth and the (n+ l)th col- 
lisions, t n , is such that the magnitude of x^ +1 is unity. 



We now seek fixed points in these equations. We as- 
sume that the time between collisions is sufficiently small 
so that the t n terms in equations ^ are negligible and 
consequently 



(10) 



during the approach to the fixed point. 

We wish to find a fixed point in the components of the 
velocity in the direction of the lines of centers. Specif- 
ically, we would like to investigate how this component 
decreases in each iteration. We can define: 



-*c -*c 7 -*c -*■ 

c n+l ' u n+l — K nX n • U 



(11) 



Taking the dot product of equations (||) respectively 
into x s n and x c n1 and using equations (|Io| ) gives 



L n+1 "n+1 



L n+2 "n+2 — *^n+l "n+1 



(12) 



1 + r . 



x n+l ' u n+l C0Sf 



(13) 



Equations (11)-([L3|) thus imply the recursion satisfied by 
the scaling factor k n 



k n +i — — — 



r 1 + ?" 



■ cos V. 



(14) 



Fixed points can be found by setting k n+ i = k n in (14), 



o ,1 + r 



k z - k—^— cos 9 + r = 0. 



(15) 



As a result, the fixed point of the scaling factor has two 
possible values 



1 + r 
k± = — - — cos 9 ± 



'1 + r 



cos 9 



(16) 



Equation ( |l6| ) is one of the major results of our study. 

In every collision, the colliding particle must approach 
particle 0. Hence, x c ■ u c must be negative in every itera- 
tion. This is possible only if k is a positive real number. 
One kind of failure arises when k is complex. Then the 
real part of the dot product will change sign infinitely of- 
ten and no fixed point can be reached. Thus, for inelastic 
collapse to occur, the quantity under the square root in 
equation ([l6]) must be positive. This positivity still per- 
mits both signs of cos 9. However, if the roots are real and 
the cosine is negative, both roots will be negative. Hence 
neither is a possible solution for inelastic collapse. The 
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only remaining possibility is that inelastic collapse may 
occur under the condition on the cosine given by equa- 
tion (|j). According to that statement, when r — > 0, 9 can 
have a value between and w. On the other hand, when 
8 = 0, the well-known one-dimensional result is recovered 
||, i.e., inelastic collapse is possible for < r < 7 — 4\/3- 
Regions (a) and (b) in figure 2 are the region of r and 9 
for which we may have inelastic collapse. 

We now consider the stability of the above fixed points. 
Stability will imply that a small change in the initial con- 
ditions will leave the system in a collapse state, or in 
other words, changes will still permit an infinite number 
of collisions. There are two collapse fixed points distin- 
guished by the values of the two multipliers k = k+ and 
k = k- . Subtracting k± from both sides of (p~4|) yields 



k n +i - k± 

kn &± 



(17) 



It is seen that fc_ corresponds to an unstable fixed point, 
while fc + corresponds to a stable one. Henceforth, we use 
k to denote the stable fixed point k = k + . 

Next, we investigate the time interval between succes- 
sive collisions. Assume that relative motion of the two 
colliding particles between each pair of collisions covers a 
distance which is very small in comparison to their radii. 
Then, we can think of the surface of the particles as flat. 
After the nth collision, the colliding particle moves away 
from the surface and covers a distance t n x n+1 ■ In 
the next step, this particle moves back over the same dis- 
tance, which is given by -t n +\X c n +2 ■ u n+2 an d reaches the 
surface. Thus, the ratio of times is 



<n+l 



"rt+1 
c ri+2 ' u n+2 



(18) 



This result may be simplified with the aid of equations 
©-© to give 



t n +i _ r_ 
t n k^ 



(19) 



Let d denote the shortest distance between two particles. 
The distance ratio equals the product of the time ratio 
and the velocity ratio, k, 



d n +i _ r 



(20) 



We performed numerical simulations of the collision 
process by considering three inelastic particles moving in 
two-dimension with random initial conditions. When col- 
lapse happens, we compared the ratios calculated from 
simulations with the predictions of equations (p"l|), (19), 
(pp|). We found excellent agreement (Figure 3), indicat- 
ing that the fixed points are attractive and indeed corre- 
spond to collapse 



IV. VALIDATION OF THE FLAT SURFACE 
APPROXIMATION 



As observed in the last section, making the approxima- 
tion (^) simplifies dramatically the original system and 
it can be described by a single ratio k n . It is as if the 
particles have flat surfaces, so that the effect of the tan- 
gential components of the velocities of particles 1 and 2 
can be ignored. This is true only if the time intervals are 
negligible, so that ( |l0| ) can be obtained from @. We will 
see when such a simplification is valid, and we will set a 
criterion for the range of validity of the approximation. 



In our way to equations (|12|)-([l3|), we neglected the 
terms like t n (u^) 2 in comparison to x n ■ u n , by using the 
approximation (|10|). As noted above, the former terms 
decrease as (p-)™, while the latter terms decrease as k n . 
Thus, the flat surface approximation is reliable only when 
r < k 3, so that terms proportional to t n can be safely ig- 
nored. This condition can be explicitly written as, 



r < 



■ cos 




(21) 



which can then be simplified into the form of condition 
(^) . The region of stability determined by this condition 
is region (a) in Figure 2. The maximum possible value 
of r for stable behavior is r c = 9 — 4-\/5. 

To this point, our calculation did not rely on the cir- 
cular geometry of the particles. The name "flat surface" 
suggests that when criterion (0) is satisfied, particles do 
not experience the curvature of their surfaces, and collide 
as if they are flat. This calculation is valid for arbitrary 
particle shapes when (|J) is satisfied. 

We also observe that when criterion (||) is satisfied, 
the time interval t n decreases faster than the radial com- 
ponent of the velocities. In such a situation, collapse 
happens so quickly that all other effects, external or in- 
ternal, have no essential influences to the process. One 
can further consider arbitrary interactions between the 
particles as well as arbitrary external fields, as long as all 
the interactions depends only on the relative positions. 
Since the particles' relative positions only change very 
little during the process of collapse, all the effects of the 
interactions on, say, particle 1 can be replaced by a con- 
stant total force acting on it which induces an constant 
acceleration. This acceleration has very little effect in 
the tangential direction since time interval is too small 
for it to change the tangential component of the veloc- 
ity. When the flat surface approximation is valid, the 
time interval is even too small to change the radial ve- 
locity component. We conclude that the the previously 
obtained fixed points are unchanged. 
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V. CIRCULAR GEOMETRY 

In the previous sections, the calculation were per- 
formed by neglecting the t n terms completely out of equa- 
tions (||) and the fixed points for inelastic collapse were 
found when the final state satisfied condition (||). After 
understanding the characteristic behavior of the collapse, 
we can do a more rigorous calculation to investigate how 
the system behaves outside the region satisfying (||). 

We now see that during the collapse process, the radial 
velocity components of particles 1 and 2 monotonically 
decrease till they vanish at the moment of singularity, 
while their tangential components approach limiting val- 
ues as the number of collision diverges. Hence, we take 
those tangential components as constants, and concen- 
trate on the radial components in the calculation of fixed 
points. 

For simplicity we study in detail one collision in the 
situation where particle 1 has a zero radius while parti- 
cle has a unit radius. Since the theory depends only 
upon the sum of the two radii, this case subsumes all 
others. Here we denote the initial instant with a sub- 
script i, and the instant before the collision with /. We 
drop the superscripts since only particle 1 is considered. 
We also assume the particles are very close, d <C 1, and 
—u-Xi <C u t , where u t is the magnitude of the tangential 
velocity component of particle 1. The collision time is 



(22) 



Immediately before the collision, the radial velocity com- 
ponent of particle 1 equals 




2du 2 t . 



(23) 



Equation (|22|) gives the quantity needed to complete 
the equation set In order to abbreviate this calcu- 

lation, we introduce an effective centrifugal acceleration. 
In the above calculation, if we view the situation in a 
frame rotating with an angular velocity of u t , around an 
axis perpendicular to both u and the line of centers, pass- 
ing through the center of particle 0, then particle 1 has 
zero tangential velocity, and the effect of the tangential 
velocity can be represented by a centrifugal acceleration 
ai = Uf. This substitution is justified by noticing that we 
can get exactly the same expressions as (^||) and (^3|) by 
using this acceleration. We do not need really use such 
a rotating frame. We use a± as an effective centrifugal 
acceleration to replace the tangential component of the 
velocity with the same effects. 

After these preparations, we are in a much clearer posi- 
tion. Particle 1 and particle 2 respectively have centrifu- 
gal accelerations a± and a2, which are all in the radial 
directions. Equation (|o|) is again a good approximation. 
Particle 1 is moving on a line, and so is particle 2. We 
can further drop the vector notation. In the following, we 



use u to denote the radial component of the velocity of 
particles 1 and 2 immediately before a collision, with its 
positive direction pointing towards the center of particle 
0. Consequently, a\ = —u\ t and a 2 — — n\ t . 
The equation set (||)-(||) reduces to 



n+l 
U n+1 



~ ru n + a n+ltn 
,i H — COS 6 



d n +i 


= —ru\ 






a n+l 


= < 


s 


c 




= a n 



u n+l ' tn 2 a «+ 1 ^™ 
+ -n s t 2 



(24) 



where d n is the distance between the spectator and par- 
ticle at the instant of the nth collision. Recall that t n is 
the time interval between the nth and (n + l)th collisions. 

One simple case can be fully carried through, the case 
a\ — ai = a. In this case, particle 1 and particle 2 are in 
a symmetrical situation so that the recursion relation of 
the system can be obtained after a single collision. 

As we did before, we use two nondimensional numbers 
k n and a n to describe the evolution of the system: 



l n+l 



and 



-a ■ t n 



(25) 



Starting from (|24|) , after some straightforward calcula- 
tion, we get following recursion relation: 




k n 



r cos0-^) - (l + r)cos6»^, (26) 



k 



1 + r 



k k 



-^(I±l C0S *-£) -(l + r )coB^. (27) 

Suppose the fixed point is (a, k), then from (p7|), we 
have 



-k 2 + ?4fcos6-k-r 
k + 1 ' 



(28) 



Since a > from its definition, the right hand side of the 
above equality must be nonegative. Thus, there must ex- 
ist two real solutions of k satisfying r.h.s.= 0. And we 
readily recover the condition (Q). 
Substituting ( f28| ) into (^), we have 



1+r 



cos 9 ■ k — r 



— (1 + r) cos 8 ■ k- 



-k 2 + i±^cos6» -k- 



(29) 
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Of course, fc+; which appeared before, are solutions of 
(^9|). From (2§) we learn that (0, k±) are fixed points. 
Let us look at other solutions of equation (^) which sat- 
isfy 



So 



1 



1 



1 + r , . , 
r H — cos # ) k 



-r = 0. (30) 



This equation has one and only one solution of k in the 
interval [0 : 1]. We denote the corresponding fixed point 
as (ceo,ko). Of the three relevant fixed points, (0, £;_) is 
unstable, and we concentrate on the stability condition 
for the other two fixed points. 

Suppose that (a n ,k n ) has a small deviation from the 
fixed point (a,k). then 



Sk n+ i 
5a n+ i 



An A 12 
A21 A22 



Sk n 
8a T: 



(31) 



where 



2 [ i-tl cos 8— y\r + (l + r) cos ( 



A12 = - 



A21 



(1 + r) cosO 



2k 2 



k 2 



A u 



Ao 



(1 + r) cos 9 - 2k 
2k 2 



Let us denote by A the eigenvalue of matrix A, hence we 
have 



A — 6A 







(32) 



where 



b 2 k 



(1 + r) cos6»- 



■ cos 6 



1 

2k 2 



1 3r 2 
+r(l+r) cosd- -^4, 

c ~ fc 5 ' 

For the fixed point (0, £;+), 



b =-nr( 1 + k +) and c= iir- 

f»y_|_ ft I 



(33) 



Ai = ^r and A 2 = 



(34) 



The point (0, £;+) is an attractive fixed point if and only 
if r < fc, . Hence we recover the condition of (J3). If we 
only require Ai < 1, we can recover condition (nT). 

For the fixed point (an, fco), we find that > if and 
only if the condition of (g) is satisfied, but that condition 
also decides the range of (r, 9) inside which we will have 



1 - b+c < 



(35) 



From equation ( p2| ) we know that one eigenvalue of ma- 
trix A is larger than 1, implying that the corresponding 
fixed point is unstable. We conclude that there are no 
additional stable fixed points beyond those which satisfy 
condition (||), and there are no stable collapses outside 
that range in the case a\ =0,2- 

We believe this is true for the general situation of 
o-i 7^ 0-2- When condition ^) is violated, our simula- 
tions show that even though the particles can be very 
close, they will get apart before having collided an infi- 
nite number of times. 

From the above calculation, and specifically equation 
( |34| ) , we can see the parameter space (r, 9) can be divided 
into three regions (Figure 2): 

a) When condition (||) is satisfied, both the eigenvalues 
of matrix A are smaller than unity. The fixed point is 
stable in all directions in space (a, k). It is the collapse 
region. 

b) When condition (||) is violated but the condition ([!]) 
is satisfied, the eigenvalue A2 > 1, while Ai < 1. In this 
region, particles can have any number of collisions before 
they might eventually separate. 

c) When the condition (|l]) is violated, both eigenval- 
ues of matrix A are larger than unity. Collapse does not 
occur. 

The above calculation is independent of the sources 
of the accelerations, which could be the interactions be- 
tween particles. The calculation confirms our previous 
argument that interactions are irrelevant in the process 
of collapse 1 . 

We also observe that in the above calculation, the cir- 
cular geometry of the particles is not essential. The cal- 
culation is also valid for arbitrary shape of the particles, 
with the corresponding radius of curvature replacing the 
radii used. Even though the radius of curvature is not 



x But when there is a sufficiently strong attraction, e. g. grav- 
ity, between particles so that the directions of the accelera- 
tions are reversed, we want fixed points with a < 0. Then the 
fixed point (a ,k ) becomes stable when condition (|^) is vi- 
olated. Inelastic collapse can happen in a much larger region 
of (r,9). 
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relevant in deciding the behavior of collapse process — the 
radii do not show up in the expression of the stable fixed 
point or condition (||), it does have some effects. The 
centrifugal acceleration, which obviously is important in 
deciding the probability of collapse, is inversely propor- 
tional to the radius. So chances are larger for collapse to 
happen when the colliding point is at a position on the 
surface with a larger radius of curvature. 



VI. CONCLUSION 

We demonstrated analytically the existence of inelas- 
tic collapse for three particles in all dimensions. At the 
last moment of collapse, the three particles have a cyclic 
behavior, which is characterized by a fixed point. We 
have established the range of the parameters for which 
the fixed point exists and the range for which it is stable. 
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